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Abstract: 

The generalised quasienergy states are introduced as eigenstates of the new integral of 
motion for periodically and nonperiodically kicked quantum systems. The photon distribution 
function of polymode generalised correlated light expressed in terms of multivariable Hermite 
polynomials is discussed and the relation of its properties to Schrodinger uncertainty relation 
is given. 
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1 Introduction 



The nonstationary guantum systems have specific integrals of motion which depend on time 
in Scrodinger representation (see, for example [1]). The integrals of the motion which are the 
operators linear in position and momentum and which depend on time in the Schrodinger 
representation have been constructed for the forced oscillator with time- dependent frequency 
and for the charged particle moving in time dependent magnetic field in [2] and [3]. The 
linear integrals of the motion for the nonstationary systems with Hamiltonian which is a 
generic multidimensional quadratic form in positions and momenta operators have been 
described,for example, in [4]. In all these systems the dependence of the parameters on time 
was the arbitrary one including both very slow adiabatic case and very fast instant change 
of parameters. 

The important partial case of nonstationary systems is the kicked systems on which the 
external forces act during very short time in comparison with all the periods of possible 
vibrations in the system under study. The aim of talk is to discuss properties of the time- 
dependent invariants for the kicked systems because several interesting phenomena such as 
extension of Casimir effect to nonstationary conditions or deformations of particle distri- 
bution functions may be demonstrated for the kicked systems. The Casimir effect is the 
phenomenon of attraction of two neutral plates between which is the vacuum of the elec- 
tromagnetic field due to the dependence of the vacuum energy on the distance between the 
plates. If one moves the plates the Casimir forces produce the work and due to conservation 
of the energy the vacuum state becomes the state with photons. So we have nonstationary 
Casimir effect which is the phenomenon of generation of photons from the vacuum. 

The electromagnetic field may be considered as the sum of modes and each mode being 
the oscillator with its own frequency. When the plates move these oscillators are just the 
oscillators with time dependent frequencies. So we treat the nonstationary Casimir effect as 
the effect existing for any system with time dependent parameters. In this case the frequences 
are these parameters (like distance between the plates). Thus the linear invariants found in [2] 
and [4] for one dimensional and multidimensional oscillators are in the case of nonstationary 
Casimir effect the integrals of motion for electromagnetic field with time varying boundaries. 

Quantum systems with energy spectra are the systems with stationary Hamiltonians.The 
dynamics of these systems is described by the transitions between the energy levels. The 
nonstationary quantum systems have no energy levels. But for periodical quantum systems 
the notion of quasienergy levels has been introduced in [5] and [6]. The main point of the 
quasienergy concept is to relate the quasienrgies to the eigenvalues of the Floquet operator 
which is equal to the evolution operator of a quantum system taken at a given time moment. 

In [7] the connection of Floquet operator and quasienergies with time-dependent integrals 
of motion has been found. Following this article we will discuss the relation of the Floquet 
operator to integrals of motion and introduce a new operator which is the integral of motion 
and has the same quasienergy spectrum that the Floquet operator has. Implicitly this result 
was contained in [4]. 

Thus,the quasienergy of the periodically kicked quantum systems may be shown to be the 
time-dependent integral of motion for such systems. If the system is kicked nonperiodically it 
is possible nevertherless to find out the time- dependent integrals of motion and to introduce 
the generalized quasienergy states [7]. It gives the possibility to extend the analisys of 



quantum chaos phenomenon which is usually based on the studying the quasienergy spectra 
of the kicked systems to the case of nonperiodically kicked systems. 

Another physical phenomenon related to the kicking is the creating the photons in 
squeezed states from vacuum due to the nonstationary Casimir effect mentioned above. 
We will discuss how the electromagnetic field (or pion or boson field of any kind) changes 
their statistical properties due to influence of the kicking. As an example the state with the 
Wigner function of generic Gaussian form will be considered as the final result of kicking. 



2 Generalized quasienergies 

Following [7] we will discuss the properties of quasienergies and generalized quasienergies for 
periodically and nonperiodically kicked quantum systems. 

If one has the system with hermitian Hamiltonian H(t) such that H(t + T) = H(t) the 
unitary evolution operator U (t) is defined as follows 

hM>=l/(*)|^0> (1) 

where \ip, > is a state vector of the system taken at the initial time moment. Then by 
definition the operator U(T) is called the Floquet operator and its eigenvalues have the form 

/ = exp(-iET) (2) 

where E is called the quasienergy and the corresponding eigenvector is called the quasienergy 
state vector. The spectra of quasienergy may be either discrete or continuous ones (or mixed) 
for different quantum systems. For multidimensional systems with quadratic Hamiltonians 
the quasienergy spectra have been related to real symplectic group ISp(2N, R) and found 
explicitly in [4]. 

We want to answer the following questions. Is the Floquet operator U (T) the integral of 
motion?The operator U (T) does not satisfy the relation 

dl(t)/dt + i[H(t), I(t)} = 0, (h = 1) (3) 

which defines the integral of motion 7(t).So,the Floquet operator U(T) is not the integral of 
motion for the periodical nonstationary quantum systems. But as it was pointed out in [4] 
any operator of the form 

I(t) = U(t)I(0)U- 1 (t) (4) 

satisfies the equation (3) and this operator is the integral of motion for the quantum system 
under study.Let us apply this ansatz to the case of periodical quantum systems. We introduce 
the unitary operator M(t) which has the form 

M(t) = U(t)U(T)U- 1 (t). (5) 

This operator is the integral of motion due to the construction given by the formula (4) for 
any integral of motion. The spectrum of the new invariant operator M(t) coincides with the 
spectrum of the Floquet operator U(T).We have proved that since quasienergies are defined 
as eigenvalues of the integral of motion M(t) they are conserved quantities. 



The given construction permits us to introduce new invariant labels for nonperiodical sys- 
tems,for example, with the time-dependence of the Hamiltonian corresponding to quasicrys- 
tal structure in time described by two (or more) characteristic times. For such systems the 
analog of the invariant Floquet operator (5) will be given by the relation 

M 1 (t) = U(t)U(t 1 )U(t 2 )U- 1 (t). (6) 

The eigenvalues of the operator M\(t) are the consrved quantities and they characteryse the 
nonperiodical quantum systems in the same manner as quasienergies describe the states of 
periodical quantum systems. This operator determines the generalized quasienergies and 
its obvious analog describes the geometrical phase which is the characteristic of special 
nonperiodical system for which the parameters of the Hamiltonian take their initial values 
after some time T [8]. 

The quasienergy spectrum of periodically kicked quantum systems may be connected 
with quantum chaos phenomenon (see [9], [10]). In [11] the integral of motion for delta- 
kicked nonlinear oscillator has been found to exist even in the case of chaotic behaviour. In 
[12] the symmetry group criterium for the periodically delta-kicked systems has been found to 
obtain either regular or chaotic behaviour of these systems. The criterium relates the Floquet 
operator spectrum to the conjugacy classes of the system symmetry group. The results of [12] 
may be applied to the nonperiodical systems too. The generalized quasienergy spectrum is 
determined by the conjugacy class of the same group to which belongs the integral of motion 
(6). For quadratic systems it is the same real symplectic group Sp(2N, R).ln the case of two 
characteristic times the classification of the possible either chaotic or regular regimes of the 
system under study coincides formally with the classification for the delta-kicked periodical 
quantum systems given in [12]. 

The geometrical phase is defined as the phase of eigenvalue of the evolution operator 
U (T) where T is the time moment at which the parameters of Hamiltonian take their initial 
values. Using similar arguments we can answer the same question as for the quasienergies. Is 
the geometrical phase the integral of motion of the nonstationary and nonperiodical quantum 
system?The answer is "yes" because such nonperiodical system has characteristic time T. 
Due to this the operator {7(£){7(T){7 _:L (i) is the integral of motion of the system. 

From the point of view of given analisys the nonperiodical systems with several char- 
acteristic times have to demonstrate similar physical properties that are usually considered 
as properties of purely periodical quantum systems. Thus, the types of chaotic and regular 
regimes of the kicked quantum systems have to be the same for both periodical and nonpe- 
riodiacal kicks and group classification of regular and irregular behaviour of kicked systems 
given in [12] may be easily extended to the case of nonperiodically kicked systems. 

3 Photon distributions 

Let us consider as nonstationary quantum system the photons in a resonator. We want to 
discuss the photon distribution function which is the distribution for the generalized cor- 
related state found in [13]. Below we follow this work. Due to the varying boundaries and 
time dependence of the refraction index of the media in a resonator the mixed squeezed 
state of the iV-mode light with density operator p may emerge. It is described by Wigner 



function VE(p, q) of the generic Gaussian form which contains 2N 2 + 3N real parameters. 2N 
parameters are mean values of light quadratures < p > and quadratures < q > and other 
2N 2 + N parameters are matrix elements of the real symmetric dispersion matrix m with 
four iV-dimensional block matrices three of which are 

mn = Op, (7) 

= 0p q , (8) 

m 22 = Oq. (9) 
Below we will use the invariant parameters 

T = Trm (10) 

and 

d = det m. (11) 
Also we will use the other invariant coefficients of the polynomial 

P(x) = det(m - xl) (12) 

where 1 is 2iV-dimensional identity matrix. 
We will introduce the notations 

Q = (p,q) (is) 

where 2iV-dimensional vector Q consists of iV components of light quadratures p±, ...,pn and 
iV components of quadratures qi, qjy. The generic gaussian Wigner function has the form 
(see, for example, [1]) 

W(p, q) = <H exp[-(2)- 1 [(Q- < Q >)m- 1 (Q- < Q >)]. (14) 

The parameters < p > and < q > are given by the formulae 

<p>=Trpp, (15) 
< q >= Tr pq, (16) 

where the operators p and q are the quadrature components of photon creation at and the 
annihilation a operators 

Due to the physical meaning of the dispersions the diagonnal elements of the matrix m must 
be nonnegative numbers, so the invariant parameter T Eq.(10) is a positive number. Also the 
determinant d Eq.(ll) of the dispersion matrix must be positive. In fact the matrix m must 
be positive-definite. 



To obtain the photon distribution function we have to calculate the probability P n to 
have n photons in the state with the density operator p. Here the vector n has N components 
rii which are nonnegative integers. This probability is given by the formula 

P n = Tr p\nxn\,ni = 0,1,2,... i = l,2,...,N (19) 

where the number states |n > are the eigenstates of the number operator afa 

afa|n>=n|n>. (20) 

The function P n may be obtained if one calculates the generating function for the matrix 
elements p mn of the density operator p in the Fock basis. This generating function is the 
matrix element of the density operator in the coherent state basis 

< p\p\a >= exp(^- - 1^) ]T f p mn (21) 
1 1 m,n=o (m!n!)2 

Here and below we use not at ions: a and (3 are iV-dimensional vectors with complex compo- 
nents and 

rx\ = n 1 \n 2 \...n N \, (22) 
a n = a^a^...a n N N (23) 

and 

oo oo oo oo oo 

E = E - E E - E • (24) 

m,n=0 mi=0 mjv=0ni=0 rijv=0 

We have 

Pn = Pnn- (25) 

The iV-mode coherent state \a > is the normalized eigenstate of the annihilation operator 

a|ct >= a\a > . (26) 

In terms of Wigner function the density operator in the coherent state representation has 
the form of 2iV-dimensional overlap integral [1] 

< P\p\a >= J W(p, q)W a(3 {p, qjdpdq, (27) 

where the function W a/3 (p,q) is the Wigner function of the operator \a >< (3\.It has the 
form [1] 

W af3 (p, q) = 2^ exp[-^ - 1^ - a(3* - p 2 - q 2 + V2a(q - ip) + v^/3*(q + ip)]. (28) 
Let us introduce the 2iV-dimensional vector 

7 = (/?*, «) (29) 



which is composed from two iV-dimensional vectors. Let us introduce also 2N- dimensional 
unitary matrix u which has four iV-dimensional block matrices 

% 

u\\ = -U12 = -iu 2 i = -iu 2 2 = ( 30 ) 

and 1 is iV-dimensional identity matrix. 

Since the integral given by Eq.(27) is the Gaussian one it may be easily calculated. So, we 
have 

< p\p\a >= P exp(-M-) exp[-i 7 P 7 + 7 Py] (31) 

where the symmetric 2iV-dimensinal matrix R has the matrix elements expressed in terms 
of the dispersion matrix m as follows 

R = -2u~ 1 {— - + l)- 1 u* + H X (32) 

The matrix T. x is 2iV-dimensional analog of Pauli matrix. The 2iV-dimensional vector y is 
given by the relation 

y = 2u*(l - 2m)" 1 < Q > (33) 
where the matrix u l is transposed one and the factor P has the form 

P = [det(m + exp[- < Q > (2m + l)" 1 < Q >] (34) 

We obtain for the photon distribution function P n the expression 



n! 



(35) 



Here the matrix R determining the Hermite polynomial is given by the formulae (32) and 
argument of Hermite polynomial is given by the expression (33). The expression (35) is the 
partial case of the matrix elements of the density operator in Fock states basis obtained in 
[14] by canonical transform method. 

The behaviour of the distributions may be wavy function of photon numbers as it is in 
one-mode case [15], [16]. 

Mean values of photons in each mode have the form [13] 

< rij >= ^(a PjPj + a qjqj - 1) + ^(< Pj > 2 + < <g). (36) 

Photon number variances are given by expression 

= \(T* - 2d, - < Q 3 > m 3 < Q 3 > (37) 

where Tj and dj are the trace and the determinant of the photon quadrature 2x2-dispersion 
matrix m 3 of the j-th mode only and the 2- vector Qj has components pj, q 3 . The correlations 
of photon numbers in different modes may be expressed analogously. 



4 Uncertainty relation and distributions 



Following [17] for one-mode case we consider the photon distribution function (35) expressed 
in terms of the series of the products of two Hermite polynomials. 

There exists the expression of the Hermite polynomial of two variables with equal indicies 
in terms of the products of two Hermite polynomials of one variable[14]. Using this formula 
one obtains from the general expression (35) the following photon distribution function 

P. = i>! £(^r<- - t)|-»H-'lg^( Jt »' ft ±g» w )|» (38) 

Here the factor P is given by the formula (34), the matrix elements of the matrix R are given 
by the formulae (32) and the numbers y are given by the formula (33). If the state is the pure 
squeezed and correlated state the determinant of the quadrature dispersion matrix is equal 
to | and the matrix element Ru = due to the formula (32). In this case only the first term 
in the series is not equal to zero and this term gives the known expression for the photon dis- 
tribution function of the pure and correlated state obtained by another method in [16]. The 
above expression is convinient to discuss the connection of the uncertainty relations with the 
quantum distribution functions. In fact the formula for the Gaussian Wigner function (14) 
formally coincides with the classical Gaussian density of probability in the particle phase 
space. Then we have to answer the question in what aspects these formulae must be consid- 
ered as essentially different ones. By intuition it is clear that the uncertainty relation must 
distinguish the classical and quantum Gaussian densities of probabilities having absolutely 
the same form (14). And in fact in addition to the usual restriction for the dispersion matrix 
for the Gaussian classical distribution function in the phase space which is the condition of 
nonnegativity of this dispersion matrix the quantum mechanics demands the inequality for 
the determinant of this matrix d > 1/4. This inequality is the Schrodinger uncertainty rela- 
tion. The formula (38) permits to relate this inequality with physically obvious condition that 
probability to find n photons must be nonnegative because the determinant of the dispersion 
matrix is the parameter on which depends the photon distribution function. We see that all 
the terms in the expression (38) are obviously positive ones except the term containing the 
number —2R\2 which may change the sign for odd powers k if it is not positive itself. It means 
that for natural condition of positiveness of photon distribution function it is nesessary to 
have inequality 

R12 < (39) 
But this inequlity is equivalent to the mentioned above inequality for the dispersion matrix 
determinant d which implies the uncertainty relation. Thus existence of the connection of 
the Wigner function (14) with the photon distribution funtion (38) is consistent only if the 
uncertainty relation holds. So we clarify the mechanism how uncertainty relation in phase 
space of electromagnetic field oscillator influences the form of photon distribution function. 

As we discussed the nonstationary Casimir effect produces the changes in statistical 
properties of the photons. For example, it changes Poisson distribution of photons to become 
the discussed distribution expressed in terms of multivariable Hermite polynomials. Also it 
deforms Planck distribution to become [14] 

n = ^- + H 2 coth^+|5| 2 , (40) 



where the correction to usual Planck distribution term contains parameters \v\ 2 and | S | 2 . These 
parameters depend on the characteristics of kicking [1], [14]. So, a generic kicking produces 
from ground state the squeezed mixed state with deformed Planck distribution and the 
number state distribution function which is expressed expliciltly in terms of Hermite poly- 
nomials of several variables. 

It is interesting to note that if to take into account a possible nonlinearity of classical 
electrodynamics we can have another reason for deforming photon distribution formulae. In 
[18] it was suggested that the possible nonlinearity of the electromagnetic field vibrations 
which must exist for very large amplitudes corresponding to high densities of the field energy 
may be considered as the nonlinearity described by the q-oscillator. It was seen that it is 
subject to nonlinear vibrations with a special kind of dependence of the frequency on the 
amplitude and the influence of such nonlinearity on Bose distribution function was evaluated. 
So, for the thermal state we have the deformed Planck distribution formula [18] of the form 

(^-T— -K 2 Tf M 1U • ( 41 ) 

e*;T — l Ki (^ e feT — lj 4 

Here the first term is the usual Planck distribution formula and its correction is proportional 
to the square of the nonlinearity parameter k. 

It means that black body radiation formula changes due to the nonlinearity of the electromag- 
netic field vibrations. So nonstationary Casimir effect and q-nonlinearity of electromagnetic 
field vibrations produce deformations of Planck distribution formula. But the temperature de- 
pendence of corrections to Planck formula is different and it gives a possibility to distinguish 
the influence of these effects. 
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